Implicit regularization beyond one loop order: gauge field theories 
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We extend a constrained version of Implicit Regularization (CIR) beyond one loop order for gauge 
field theories. In this framework, the ultraviolet content of the model is displayed in terms of 
momentum loop integrals order by order in perturbation theory for any Feynman diagram, while 
the Ward-Slavnov- Taylor identities are controlled by finite surface terms. To illustrate, we apply 
CIR to massless abelian Gauge Field Theories (scalar and spinorial QED) to two loop order and 
calculate the two- loop beta- function of the spinorial QED. 
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I. INTRODUCTION 



Alternatives to dimensional regularization (DR) 
need to be constructed for quantum field theoret- 
ical models which are well-defined only in their 
physical dimension. For instance, in supersym- 
metric theories the analytical continuation of the 
space-time dimension from 4 to d for the vector 
fields leads to a mismatch between the fermionic 
and bosonic degrees of freedom. This gives rise 
to a breaking of the supersymmetric relations. A 
fully consistent and symmetry preserving regular- 
ization and renormalization framework, valid to ar- 
bitrary loop order, has not yet been constructed. 
Instead, a more pragmatical approach is employed 
by using Dimensional Reduction (DRed) 1]. Al- 
though DRcd is mathematically inconsistent and 
cannot be extended to arbitrary loop order, it can 
be successfully applied if it is judiciously used in or- 
der to preserve the supersymmetric Slavnov- Taylor 
identities of the model. In principle, an invari- 
ant regularizaton scheme could be precluded by 
adopting an strategy in which symmetry breakings 
are compensated by symmetry restoring countert- 
erms, once genuine (physical) quantum symmetry 
breakings have been discarded in the underlying 
renormalizable model. In practice, a more efficient 
approach would be provided by a manifestly su- 
persymmetric and gauge invariant regularization. 
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From a more phenomenological perspective it is 
expected that the Minimal Supersymmetric Stan- 
dard Model (MSSM) can be probed at LHC with 
an accuracy at the percent level in measurements 
of electroweak precision observables 0. This, on 
the theory side, is realized through virtual effects 
of supersymmetric particles. The theoretical eval- 
uation of such observables must be performed at 
least to two loop order both in the SM and MSSM 
so that the evidence of new physics becomes un- 
raveled. Clearly, an invariant regularization and 
renormalization procedure is important, both to 
eliminate inconsistencies and to reduce the number 
of adjustable parameters in the calculation, partic- 
ularly when a supersymmetry breaking mechanism 
is taken into account in more phenomenological 
physical models. 

Implicit regularization (IR) technique 0)-[I3 
is a relatively new momentum space framework 
which operates in the physical dimension of the 
theory. Under the assumption that some regulator 
is implicitly being used in a divergent amplitude, 
it is possible to algebraically separate the diver- 
gences in terms of momentum loop integrals which 
are independent of external momenta, and may be 
directly absorbed in the definition of the renorma- 
lization constants. In this framework, a minimal 
subtraction consists in defining the renormaliza- 
tion constants using mass independent loop inte- 
grals. This is achieved by exchanging the mass 
dependence for an arbitrary non-vanishing scale 
by means of a regularization independent identity. 
Such scale parametrizes the freedom of separat- 
ing the divergent part of an amplitude and plays 
the role of renormalization group scale. IR can 
be generalized to arbitrary loop order 0] , [H[ . As- 
suming that the theory has been renormalized to 
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n — 1 loop order, it is always possible to cast the 
ultraviolet behavior in n loop order as a momen- 
tum loop integral. In other words IR is compatible 
with the BPHZ forest formula which defines the 
set of subtractions to remove subdivergences. The 
derivatives of the renormalization constants which 
are useful to calculate renormalization group func- 
tions can be displayed as loop integrals as well. 
Because the loop basic divergent integrals are not 
evaluated, nor is the integrand of the amplitude 
modified, such procedure neatly separates the reg- 
ularization dependent part (loop momentum inte- 
grals) from the finite part (which can be evaluated 
by standard methods [181 ]) in a unitary preserv- 
ing fashion. IR has been shown to be a symmetry 
preserving framework. We have verified that sym- 
metry breaking terms stem from well defined finite 
differences between loop divergent integrals of the 
same superficial degree of divergence. Those dif- 
ferences are arbitrary (regularization dependent) 
numbers which can be written as surface terms and 
appear matter-of-factly in the process of defining a 
complete set of basic divergent integrals within IR. 
Yet such symmetry breaking surface terms can be 
handled by either choosing finite counterterms to 
restore the symmetry of the renormalized Green's 
functions, or by adding local counterterms into 
the Lagrangian so that the Slavnov- Taylor iden- 
tities are confirmed, an invariant scheme is cer- 
tainly more appealing. A constrained version of 
IR (CIR) in which those surface terms are set to 
zero ab initio has been shown to automatically pre- 
serve abelian and non-abelian gauge symmetry to 
one loop order and supersymmetry to three loop 
order 0,0,0, Ell, Tt is noteworthy that 
such surface terms vanish if they are computed in 
Dimensional Regularization which is notoriously a 
gauge invariant scheme though it breaks supersym- 
metry due to the space-time dimensional contin- 
uation. CIR was mapped onto Constrained Dif- 
ferential Renormalization [l9j], a coordinate space 
perturbative renormalization program which auto- 
matically delivers Green's functions as functions 
of a single mass (renormalization group) scale 
that fulfill the corresponding Ward-Slavnov- Taylor 
identities of the underlying model. Although Dif- 
ferential Renormalization has been successfully ap- 
plied to higher order calculations, its constrained 
version has not been established yet [20| . 

It is instructional to discuss the role of the sur- 
face terms that appear within IR in respecting the 
symmetry of a model on the quantum level. We 
have verified that setting surface terms to zero 
automatically ensures momentum routing invari- 



ance in a Feynman diagram 0,0. In other words, 
shifts in the integration variables will be allowed. 
On the other hand it is well known that the di- 
agrammatic proof of the Ward-Takahashi identity 
for correlation functions in QED involves a shift in 
the integration variable [2l| which is always legit- 
imated in Dimensional Regularization. Moreover 
when anomalies (quantum symmetry breakings) 
are present they manifest themselves as a break- 
ing of momentum routing invariance in perturba- 
tion theory [12] ■ IR correctly displays the anomaly 
among the Ward identities in a democratic way 
just as a sound regularization framework should 
by letting physics to fix the identity to be vio- 
lated 0,[1|,E3|- A good example of this feature 
is the AW triangle anomaly: the conservation of 
the vector or axial-vector current depends on the 
context it arises as the answer is not intrinsic to 
the AVV triangle graph. In such cases, as well as 
when radiative corrections appear to be finite and 
undetermined 23], surface terms are left as free 
parameters to be fixed on physical grounds. 

The idea of associating surface terms and vari- 
able of integration ambiguities to the preserva- 
tion of conventional and super Ward-Takahashi- 
Slavnov Taylor identities is not new. In a pro- 
gram called Pre- Regularization [24] , it was verified 
in some examples that a particular choice of the 
momentum shift parameters in conjunction with 
dimensional regularization could be used to pre- 
serve gauge and supersymmetry. IR goes some- 
what in the opposite direction. The vanishing of 
surface terms which also leads to momentum rout- 
ing invariance turns up to preserve gauge and su- 
persymmetry systematically. 

In a previous contribution, we have shown that 
IR can be applied to arbitrary loop order by mak- 
ing use of the BPHZ forest formula to remove sub- 
divergences [H. Basic divergent integrals, as well 
as new surface terms are consistently defined to 
n th loop order. Given the discussion above, it is 
natural to ask whether CIR respects gauge sym- 
metry to arbitrary loop order. Bearing in mind 
the Ward-Takahashi identity which states that if 
M. = e^M^ is an amplitude contributing to the 
S-matrix (e M being the polarization vector of an 
external gauge boson), then k^M^ = 0, for this ex- 
ternal gauge boson with momentum k^, the answer 
to this question is positive. Such statement stems 
from the fact that gauge invariance and unitarity 
require that the non-transverse degree of freedom 
of any massless gauge boson should not contribute 
to physically measurable quantities. The general 
diagrammatic proof of the Ward-Takahashi iden- 
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tity assumes the possibility of a shift in the inte- 
gration variable. Within CIR this is achieved auto- 
matically by setting surface terms to zero from the 
start. In other words, an invariant regularization 
and renormalization scheme emerges by systemat- 
ically setting the surface terms to zero. To illus- 
trate, we use scalar and spinorial QED as testing 
grounds. 

The paper is organized as follows: In section 
II, we present the procedure of CIR to n-loop or- 
der, discussing the strategy to subtract the diver- 
gent contributions. We also illustrate the proce- 
dure with a three-loop order calculation in the con- 
text of spinorial QED. Section III is reserved for a 
discussion on the preservation of gauge symmetry 
and unitariry by Constrained Implicit Regulariza- 
tion. In Section IV, we apply CIR to the calculus 
of the vacuum polarization tensor in scalar QED. 
Section V is dedicated to spinorial QED, also to 
two- loop order: the calculation of the vacuum po- 
larization tensor, of the electron self-energy and of 
the electron-photon-electron vertex correction are 
performed. The Ward-Takahashi identities are ver- 
ified and the beta-function is calculated. Section 
VI is dedicated to the conclusions. Some technical 
points are left to the appendix. 



II. N-LOOP CONSTRAINED IMPLICIT 
REGULARIZATION 

In this section we state the basic steps of CIR 
to n-loop order. The goal is to identify the typical 
divergence of the n th order and the finite part of an 
amplitude once the model has been renormalized 
to (n — l) th order. 

1. In order to give mathematical rigor to any 
algebraic manipulation performed in the am- 
plitude, we implicitly assume that a regular- 
ization has been applied. It can be main- 
tained implicit, the only requirement being 
that the integrand of the amplitude nor the 
dimension of the space-time is modified. A 
good one is the simple ultraviolet momentum 
cutoff. As we will see, the possible symmetry 
violations are restored by construction. Af- 
ter performing symmetry group operations, 
we cast the momentum-space amplitude as a 
combination of basic integrals. 

2. Once the basic integrals are obtained, the di- 
vergent part is separated as loop divergent 
integrals, which are obtained by applying re- 



cursively the algebraic identity, 

1 _ 1 

(p — k) 2 — to 2 (k 2 — m 2 ) 

P 2 - 2p ■ k 

(k 2 — m 2 ) [{p — k) 2 — to 2 ] ' 

in all propagators that depends on external 
momenta, until the divergent part is freed 
from the external momentum dependence in 
the denominator. This will assure local coun- 
terterms. 

3. Loop divergent integrals for which the inter- 
nal momenta carry Lorentz indices are ex- 
pressed as functions of basic divergent in- 
tegrals and surface terms (the n-volume in- 
tegral of a total n-divergence) . Such sur- 
face terms are regularization dependent (e.g. 
they vanish in dimensional regularization) 
and are connected to momentum routing in- 
variance: local momentum routing depen- 
dent terms show up in the evaluation of an 
amplitude multiplied by (dimensionless) sur- 
face terms, which may break gauge symme- 
try. CRI assumes that such surface terms are 
canceled by local symmetry restoring coun- 
terterms. In practice, this is automatically 
realized by setting them to zero from the 
start. Care must be exercised when quan- 
tum symmetry breaking occurs. In this case 
the surface terms should be considered as 
finite arbitrary parameters to be fixed on 
physical grounds. That is because anoma- 
lies are somewhat connected to momentum 
routing dependence in Feynman diagram cal- 
culations [221 ]. 

4. The basic divergent integrals which encode 
the ultraviolet behavior of the amplitude 
need not be evaluated. We adopt the fol- 
lowing notation for basic divergent integrals 
of n th loop order: 

(2) 

r(») f 2\ f 1 ryOn— 1)/,2 2 \2\ 

I lU m )=j k Jj^r/o '(*,m,A) 

(3) 

for logarithmically and quadratically basic 
divergent integrals, etc. . Hereafter the su- 
perscript A indicates that integral is regular- 
ized and J, = J d 4 k / (2tt) 4 . The functions 
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Zq™ ^ are typical finite terms of the (n — 
l) 4,l -loop order. A is a non-vanishing arbi- 
trary parameter with mass dimension which 
appears when we define a mass independent 
scheme as we explain below. At one-loop 
order the basic divergent integrals in four 
dimensional space-time reads Ii g(m 2 ) = 
J k (k 2 - rn 2 )- 2 , I quad {m 2 ) = J k (k 2 - m 2 )' 1 , 
etc. 

5. The subtraction of the subdivergences is per- 
formed in the light of the BPHZ Forest For- 
mula [25J . We adopt the version of the for- 
est formula which is analogous to the ordi- 
nary counterterm method. The subtraction 
operators are translated into local countert- 
erms which substitutes the subgraph. This is 
implemented via a recursion equation which 
involves disjoint renormalization parts only, 
described in textbooks [26j]. The basic di- 
vergent integrals can be subtracted as they 
stand in the definition of renormalization 
constants. A minimal, mass independent 
scheme is defined by substituting m 2 with 
A 2 7^ using a kind of scale relation which is 
regularization independent. This is achieved 
by means of the identity, 



1 



1 



{k 2 ~m 2 ) w (k 2 -\ 2 Y 



(A 2 - m 2 ) 



E 



^ {k 2 -m 2 Y{k 2 - A 2 )"-^ 1 



(4) 



For a basic logarithmic integral in one loop 
order, such scale relation is 



h og {m 2 ) = h og (\ 2 ) + b\n 



A 2 



(5) 



where b = 



(4tt) 



The object to be subtracted 



is Iiog(^ 2 ), A playing the role of renormaliza- 
tion group scale. For infrared safe massless 
models a systematic cancelation of ln(m 2 ) 
coming from ([5]) and from ultraviolet finite 
part will render the amplitude well defined 
as m 2 — > 0. 

6. The remaining ultraviolet finite integrals are 
evaluated as usual using Feynman parame- 
ters or other methods in momentum space 

m- 

We start off with some examples. In order to 
illustrate the role of surface terms in the preserva- 
tion of symmetries, we recall the one-loop massless 



QED vacuum polarization tensor, as calculated in 
[5|, with arbitrary momentum in the internal lines, 
k\ and k 2 . It is written as: 

n M „ = n(p 2 )(p^p^ -p 2 g^) 
1 



+ 4 I a-ig^ --(k 1 + k 2 )a 2 g^ 

+ i(fcf fcf + k$ki + k<?k%)azg {lxv g a(}} 
- (fci + k 2 ) a {ki + k 2 )^a 2 g va 



-(fcffcf + k%k%)g^a 2 g af} 



(6) 



In the equation above, p = k\ — k 2 is the external 
momentum and 



7 iog (A 2 )-61n( -^j+-6 



2 

V \ 



with b = z/(47r) 2 , includes the basic divergent inte- 
gral (in which we have used[5]). We have chosen the 
massless limit just for the sake of simplicity. Now, 
the momentum routing dependent terms, which 
lead to violation of gauge symmetry, are propor- 
tional to the aj's, namely 



otig^v 



a 2 g^v 



9n» 



k (k 2 — m 2 ) 2 



(7) 



9^ 



_ 4 [ , (8) 



k (k 2 — m 2 ) 2 J k (k 2 — m 2 ) 3 



and 



a i9{nv9ccP} = 9{pLv9a()} 



24 



k (k 2 — m 2 ) 2 
kfikyk a kp 



k (k 2 — m 2 ) 4 



(9) 



These parameters are surface terms. It can be eas- 
ily shown that 



&29/1V 



atig^v 



A d ( K 



k dk^ \ (k 2 — m 2 Y 



d ( k v 



k dk^ \ (k 2 — m 2 ) 



(10) 



(11) 



and 



r A d 


Ak^ki/ka 


L Ok? 


(k 2 — m 2 ) 3 



9{nu9a0}( a 3 - OL 2 ). 



(12) 
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FIG. 1: a contribution for the three loop massless QED 
self-energy 



As a second example illustrating the steps of CIR 
at n-loop order, let us now evaluate the nested 
three-loop contribution to the fermion propagator 
in massless QED as depicted in fig. [TJ 

We begin with the one-loop contribution, that 
reads 



iX(p) = -e 2 



(13) 



Ik k 2 {p-k) 2 
and which may as well be written like 

i£(p) = 2e 2 j a I a , (14) 
where I a is a linearly divergent integral 



I a = 



k k 2 {p-kf 



(15) 



In order to calculate this basic divergent integral, 
we apply the identity ([T]) twice to obtain 



I a = 



k a 



1 



k (k 2 — to 2 ) [ (fc 2 — m 2 ) 
p 2 - 2p ■ k (p 2 - 2p ■ k) 2 



(k 2 — to 2 ) 2 (k 2 — m 2 ) 2 [(j> — k) 2 — to 2 



(16) 



in which the limit to 
of the calculation. Thus 



will be taken in the end 



4 = 



k k 



k (fc 2 — TO 2 ) 3 

k a (p 2 -2p-k) 2 



k (k 2 — to 2 ) 3 [(p — k) 2 — to 2 



I a being a finite integral. In the limit to 2 — > we 
obtain, taking ([5]) into account 



In 



Ilog(^ 2 



bill 



El 

A 2 



2b 



y "257,0 



(18) 



where a 2 is the surface term defined in ([5]) that 
will be set to zero. Thus 



2b\ (19) 



/-(/-) = e 2 ^/ io9 (A 2 )-&ln(-^ 



Now consider the nested two-loop self-energy. Ac- 
cording to the BPHZ forest formula the subtrac- 
tion of the one loop subdivergence amounts to re- 
place the inner diagram with its finite part, namely 



iEf\p) = -ibe 4 [ 

Jk 

= 2ie 4 bj a 



k k 2 { P -k) 2 



In 



k 2 , 

A 2 2 



21 -I (2) 



(20) 



In the equation above 



/(2) = 



1 1-^ 

k k 2 ( P -k) 2 \ A 2 



Again we apply identity (|TJ) twice, to obtain 



j(2) 



In 



(21) 



1 



k (k 2 — m 2 ) \ A 2 / [ (k 2 — to 2 ) 
p 2 - 2p ■ k (p 2 - 2p ■ k) 2 

(k 2 — m 2 ) 2 (k 2 — TO 2 ) 2 [(p — k) 2 — to 2 



(22) 



Following the same steps as before, we have: 



Ii 2) = W 



k k 
™fA"'a 



k (k 2 — TO 2 ) 3 

k a (p 2 -2p-k) 2 



In - 



k (k 2 — m 2 ) 3 [(p — k) 2 — m? 



In 



k 2 



A 2 







k dk^ \ (k 2 — to 2 ) 2 



In - 



9/J.O 



1 



k (k 2 — to 2 )' 



k (k 2 — to 2 ) 2 



In - 



k 2 - 



(23) 



' 1 j _|_ i a (j) 2 m 2 ) (17) la being a finite integral. A cancelation of diver- 

k (k 2 — to 2 ) 2 J gent logarithms between the ultraviolet finite and 



G 



divergent parts as m 2 — > takes place to render 
the amplitude well denned, namely 

r(2) _ Pa / r(2),],2\ , 1 j ( X 2s 



In 2 f -I 



-In - 



A 2 



-3 



(24) 



in which we have adopted the following notation 
for the basic logarithmically divergent two-loop in- 
tegral for massless theories, 



4> 2 ) 



i 



k [k 2 — m 2 ) 2 



In 



and we have used the two-loop scale relation 



(25) 



Jg(m a )=-0* a ) 



In 2 



A 2 



In 



A 2 



(26) 



Altogether we get, for the two loop nested electron 
self energy 



i^\p) = l e%M-P i :>(\ 2 ) + ^I log (\ 2 ) 



log V' v J n 2 lo 9 y 



In 2 



A 2 



-bin 



?1 
A 2 



b .(27) 



Finally, we calculate the three-loop nested self- 
energy. Again, according to the forest formula, to 
compute the divergence of the order, the subtrac- 
tion of the subdivergences amounts to plugging the 
finite part of the two-loop subgraph into the total 
three-loop amplitude, i.e. 



iE< 3) (p) 



-le 



k 4 (p - kf 



(28) 



where the tilde stands for the finite part of (|2"7|) . 
Therefore we are able to write 

l ^\p) = -e% 2 r(l { a 3) -5I^+5I a ), (29) 



with 



1^ 



k k 2 (p-k)'- 



In 2 



A 2 



(30) 



We proceed as before. In the integral above, we 
isolate the divergence with the help of identity (fl]) , 



which is applied twice to get 



I { a 3) = 2pf 



In 2 



k (k 2 — to 2 ) 3 \ A 2 
k a (p 2 -2p-k) 2 



k {k 2 — m 2 ) 3 [(p — k) 2 — m? 



A 2 



A 



( k a \ n 2 (- k2 ~ m2 



2 1 J k dk^ y [k 2 — m 2 y 2 \ A 2 

1 

r 



.j(3) 



where 



' ~ > k (k 2 -m 2 ) 2 



In' - 



7 2 2 

k — m 



A 2 



(32) 

The complete result of (f3TJ)) is achieved when the 
surface term is set to zero, the three-loop scale re- 
lation, 



2 * 

m 



In" — +21n — 



a 2 ; 1 ~^va 2 ' ' (33) 

is used and the finite part is calculated. We have 

^ = f {4 3 i(A 2 )+4 2 ](A 2 ) + iwA 2 ) 



\ , b m 2 /" P 



3 ln V A 2 / ' 2" V A 2 



(34) 



Finally, the three-loop contribution to the self- 
energy of fig. [T] is given by 

^ 3) (p) = -^{/S(A 2 )-<(A 2 ) 



+3/ ;os (A 2 ) - £ In 3 ( ) + 36 In 2 ( - 



2 / P_ 

A 2 



(35) 



The procedure we described in this section will be 
used henceforth to compute the two-loop scalar 
QED and spinorial QED. As it will be shown, the 
elimination of surface terms plays essential role in 
satisfying the Ward identities. 
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III. DISCUSSION ON GAUGE 
INVARIANCE AND UNITARITY 



In this section, we show that CIR is built so as to 
preserve gauge symmetry and, as a consequence, 
unitarity. We are concerned here with QED, al- 
though this argument can be easily extended to the 
non-abelian case. As discussed in the previous sec- 
tion, the prescription to be followed assumes that 
some regularization has been applied. As we will 
see, this regularization does not necessarily have to 
preserve gauge symmetry, since gauge invariance 
is restored by construction. So, we assume here 
that the regularization is a simple ultraviolet cut- 
off, which does not modify neither the dimension 
of space-time nor the integrand. The procedure 
we are adopting is the traditional diagrammatic 
proof of gauge invariance of the text books (see, 
for example, ref. [21j). with the addition of the 
regularization. 

Let us consider the amplitude, M.(k) = 
e tl {k)M.^(k), for some QED process, involving an 
external photon with momentum k. Then the 
Ward-Takahashi identity requires that 



so that 



k^M^k) 



0. 



(36) 



Since the external momenta are not necessarily on- 
shell, some of the contributions to the r.h.s of (|36[) 
are non-null , but they do not contribute when the 
5-matrix element is extracted. We can begin by 
any specific diagram that contributes to M. If 
the external photon 7(fc) is removed, we get a new 
diagram that contributes to a simpler amplitude, 
M.q. If now we sum over all possible insertions of 
this photon in Mo, we must obtain eq. (|36l) . 

Concerning the possibilities of insertion, we have 
two situations: the insertion can be made in an 
open fermion line or in a closed fermion line. For 
the case of an open fermion line, there is nothing 
different from the procedure of text books. These 
contributions are the ones which are non-null in the 
r.h.s of (|36p . but which do not contribute when the 
S'-matrix is considered. So, we dedicate ourselves 
to the case of insertions in a closed fermion line. 

Let us consider this line has n vertices. The 
fermion momenta for this line are p\, P2, • ■ ■ , p n - 
If the insertion is done in a point with momentum 
Px, all the fermion propagators after the insertion 
have the momentum increased by k. We use the 
identity 



t = [{j/>x +#-m)-tf 1 - m)] 



(37) 



^ x + ^ — m j) x — m 

i i 



j> x — m ^ x + ^ — : 



(38) 



The regularized integral over the loop momentum, 
Px, takes the form 





■'N 
■ f A 




% \ tr 

Jpi 


"( 


i 


i> x -m 



i 



r 



]& 2 + ft — m 



r 



r 



(39) 



When the insertion is performed in all other points 
and the contributions are summed, there occurs 
many cancelations and only two terms survive. We 
are left with 



i-A 

il tr 
' pi 



I'n 



i 



r 



v/'l 



h- 



tx + t- 



-r 



(40) 



The two terms differs by a shift p\ — > p\ + k. If 
the integral is divergent, the cutoff regularization 
does not allow a shift without compensating with 
a surface term. So, after shifting we get 



rp/J.1 



(41) 



with S^ 1 ' being the surface term that causes 
violation of the Ward-Takahashi identity. The pre- 
scription of Constrained Implicit Regularization 
is (whenever there are no anomalies) the restora- 
tion of the symmetry by eliminating all the surface 
terms. CIR provides a simple algorithm to identify 
all the surface terms. The procedure is explained 
by the rules 2 and 3 of the previous section: iden- 
tity (TT]) is recursively used in all propagators which 
depends on the external momenta until the diver- 
gent part can be written in terms of loop integrals. 
The surface terms come from some of these loop 
integrals which have Lorentz indices. 

Concerning the unitarity, we call the reader's at- 
tention to the fact that, since the Ward-Takahashi 
identities are respected, the amplitude under con- 
sideration, after the elimination of the surface 
terms and subdivergences, can be written as 



At. 



L, 



(A A + A) 



(42) 

In the equation above, the tensor L^...^ has the 
necessary structure so as to preserve gauge invari- 
ance, A A is the divergent part (cutoff dependent) 
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forest formula, we can write 



FIG. 2: two loop polarization vacuum tensor correc- 
tions in scalar QED 



and A is the finite part. The divergent piece con- 
tains only loop integrals, which will ask for local 
counterterms. From (|42|) it is clear that the coun- 
terterms are also symmetric. So, unitarity is pre- 
served. 



IV. EXPLICIT CALCULATIONS: WARD 
IDENTITIES IN HIGHER ORDER 
CALCULATIONS 

Ward identities are relations which need to be 
satisfied by the Green's functions in order to guar- 
antee gauge invariance, unitarity and renormaliz- 
ability. From the perturbative standpoint, such re- 
lations must be satisfied order by order and should 
not be spoiled by an invariant regularization. In 
this section, we illustrate how CIR is a gauge in- 
variant scheme by verifying some Ward identities. 



A. Vacuum polarization tensor in Scalar 
QED 

At two loop order, we have four non-vanishing 
contributions, in massless scalar QED, for the 
vacuum polarization tensor, as shown in fig. [2l 
Tadpoles diagrams do not contribute because 
quadratic divergences can be shown to vanish in 
massless theories within IR [|| just as in Dimen- 
sional Regularization. 

In the diagram with nested subdivergences, fig. 
2a, we can notice that the subdiagram is the self- 
energy correction for the scalar particle at one loop 
order. After applying the Feynman Rules for scalar 
QED and subtracting the subdivergences with the 



-2e 2 k 2 b\n 



d 4 k {p + 2k)^{p + 2k) l , 
(2tt) 4 [(p + fc) 2 ](fc 2 ) 4 
k 2 



AV 



(43) 



The amplitude can be displayed in terms of typical 
basic integrals: 

Ii;„(p 2 ) = 2ie%p^ 2 Xp 2 )-Apj( 2 Xp 2 ) + 

4/^) (p 2 ) - 2p^pJ(p 2 ) - SpJ^p 2 ) - 

(44) 



8/ F (p 2 )], 



where 



rW r(0 rW = /" A M^M^ . j_i 



A 2 



(45) 

For i = 1 we simply write /, 1^, etc. 

Next, we treat the contribution coming from fig. 
2c. The subdivergences in this case are propor- 
tional to basic quadratic divergences and, so, van- 
ish for the massless case. This amplitude reads 



nUp 2 ) 



-4«e 4 5 /iI , 



1 



k,i (p-mm-k) 



A 



1 



k (P~k) 2 
k 2 



x [ Ii og (X z ) - b\n [ -—)+2b 



(46) 



where we have used the result of the one-loop ba- 
sic integral / in the internal momentum I. In the 
equation above, taking into account that shifts are 
allowed in Implicit Regularization, in our case only 
the term in the logarithm survives. This is be- 
cause the others are proportional to quadratic di- 
vergences. So, we are left with 



1 



{P - k) : 



In 



fc 2 
A 2 



= 4 4 6eV^ (2) - (47) 
The next contribution is the one from fig. 2d: 



II* (p 2 ) = 2ie 



l f A (2fc-p) M (2fc-Q„ 

k ,k 2 (k~p) 2 l 2 (k-l) 2 - 1 1 



For this diagram, the elimination of the subdiver- 
gences is carried out by considering only the finite 
part of the integral in the internal momentum I. 
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So, using the results of the one loop integrals, we 
have 



2ibe [ 61^ - -p^PvI 



31 



(2) 



(2) 



(49) 



In the expression above it was also taken into ac- 
count that 1^ — p^I /2. 

For the overlapped diagram of fig. 2b, we use 
the forest formula to subtract the subdivergences. 
So, we have 



it, 



IT 



61 



2LT 



bCT 



l±V fill ' fJ,U 

For the calculation of the first term, we have 



(50) 



n 



61 



-ie 



(2fc-p)„(2fc-p)„M 
kl k 2 l 2 (k-l) 2 (p-l) 2 {p- kf 



with 



(51) 



M = (k + l) ■ (k + l - 2p) 

= 2{(l-pf + l 2 -(k-lf -p 2 } . (52) 

Thus, it is possible to write 



n 



61 



-2 ie 4 {V(4^ 2 - W /°) 



+ [/ ;o3 (A 2 ) + 2&] [AI^-p^I] 



-2b 



2I^J-3p^+p^I^]}, (53) 



where 



1 h h h hi 



tO tO tOI t02 _ 

' J kl k 2 P(k-l) 2 ( P -l) 2 {p-k) 2 ' 

(54) 

Notice that the integrals above are symmetric in 
the exchange k <-> I. They will appear hereafter. 
Moreover 1° and I® — p^I /2 are finite integrals. 
For the counterterms, we have 



IlfJ = ie 4 I log (\ 2 ) [41^ - Vtl p v I\ , 



(55) 



and, so 



n b ,„ = -2 l e 4 {V(4^ 2 -p^/°) 

2lj l 2 J-3p^+p^I^]}. (56) 



\1V 

-2b 



We see that all the nonlocal divergent terms have 
been eliminated. Adding all the contributions, 



IF 



2H^ 



n 



■n?. 



4n 



l_i is ' 



(57) 




FIG. 3: Spinorial QED contributions for the polariza- 
tion vacuum tensor at two loop order 



(2) 



we are allowed to write 

IV = ie% {-16 Pfl I^ + S Pt ,p v I^ + Ag^ 2 

(58) 



-^[4^ 2 -W°]} 



Writing explicitly each basic integral, using the re- 
sults of the appendix, we have 



n 



- (P 9n» -P^Pv) \ 4/; os (A ) - 46 In ( ^ 2 



P 



2 tO 



,b 2 TT 2 



P 



9 



401 
"3J 



(59) 



So, we have obtained explicitly the transversal 
form of the vacuum polarization tensor in two-loop 
scalar QED, as required by gauge invariance. 



B. QED Ward Identities to two- loop order 

We now apply CIR to spinorial QED to two loop 
order. First, we shall obtain the vacuum polariza- 
tion tensor, which has only two contributions, as 
depicted in fig. [3l where we also present the cor- 
respondent counterterms, that eliminate subdiver- 
gences. 

In the first diagram, the nested subdiagram is 
the one-loop fermionic self energy correction. Tak- 
ing only the finite part of the subdiagram, we will 
correctly eliminate the subdivergences. We write 



TL= I trU-ie^^C-ie^-iSOO^ 



= l6e 4/ tr( 7 (^-^) ^ 



k 2 (p — k) 2 



k 2 
A 2 , 

(60) 



After some straightforward algebra, we have 
T« v = 4ie%[-21$ + 2p^ + 4V 
-2p,pJ +p 2 ^(2I - /< 2 >)}61) 



10 



The second contribution corresponds to an over- 
lapped diagram, for which we have to consider the 
two equal counterterms represented in fig. 3b: 



rpb 
[If 



rpbCT 



(62) 



For the first part, we have 



p61 



tT{fuh p hM-t>h P (l-t>)} 



k k 2 i 2 (k-i) 2 (p-k) 2 ( P -iy 



(63) 



which, in terms of basic integrals, gives us 

T« = -ze 4 {2 (l log (X 2 ) + 2b) (-I^+PpPuI) 
+2b(l^-2p^)-\p^(I) 2 



V^[-4(/ ;o3 (A 2 ) + 26)/+p 2 / c 
+ (/) 2 +4W( 2 ']}. 



(64) 



For the counterterms, we have 
2T% T = 8ie% 



log {\ 2 )[2I^-p^p v I + p 2g -fl). 

Now, using the results of the appendix, we add 
these contributions in order to obtain the complete 
spinorial QED vacuum polarization tensor: 



T» v = ^e 4 6( W - W 2 ){|/ 2 s (A 2 )-36/g(A 2 ) 



+ | &W A 2 )-^ln(-g) 

+ h- P 2 i°}. 



(66) 



Now, we verify the Ward Identity which relates 
the fermion-photon vertex and the fermion self- 
energy corrections. In every order of perturbation 
theory, they must obey 



(67) 



where q = p — p' is the momentum of the exter- 
nal photon. We start by evaluating the two-loop 
fermion self-energy, for which picture [5] represents 
all non-null contributions. 

The amplitude of figHh- has been already evalu- 
ated in section II and is given by 



iI>W(p)=-2ie*lrf (l {2} - P ,l) 



(68) 



We turn ourselves to the next two-loop self-energy 
contribution, which is represented by figHb, where 



(a) 



(b) 



FIG. 4: contributions for the fermionic self energy at 
two loop order 



we can notice that the subdivergence is due to the 
presence of the one-loop vacuum polarization ten- 
sor. The expression for this diagram is 



1 - L i p m l 



(69) 



where ILpJ is the finite part of the one-loop vacuum 
polarization tensor. It is simple to show that 

^i 2) (p) 

= Ue% t In 
3 k k\p-k) 2 

4 



A 2 



-b^ {2p,I^ - 2J< a > 



+ ~P»I • (70) 



The overlapped diagram in fig- St, with coun- 
terterms included is given by 



*Sl 2 ) = *El 2) + 2.E?> 



cCT 



(71) 



The contribution of the first term can be easily 
shown to be 



4»cV {-^[^o S (A 2 ) + 2b]p,I(p 2 ) + bl^(p 2 ) 

(72) 

The two counterterms give identical contributions 



2i£<g T (p) = 2ie^I log (\ 2 )I, 



(73) 



and, so 

iVfXp) = Aie^-p^p 2 ) + j( a V)l- (74) 
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(a) 



(b) 





(d) 



(e) 







(c) 





(f) 



(g) 



FIG. 5: contributions for the vertex function at spinorial QED 



Summing these three self-energy amplitudes, we 
finally obtain 



£ (2) b) = 

4 4 r A w - m 

3 J k k^(p-kf 



•» 1 4 



+ ^ 4 ^{-4 2) +4^/ (2) + ^/}. (75) 



For the evaluation of the vertex corrections, we 
would like to call the reader's attention to the fact 
that, although Implicit Regularization does not 
commute with the contraction of the amplitude 
with the metric (in this process we would loose 
control of the surface terms, which possibly are 
symmetry violating terms), it commutes with the 
contraction with an external momentum. For this 
reason, for the purpose of verifying the Ward iden- 
tities, we shall evaluate the vertex corrections con- 
tracted with the momentum of the external pho- 
ton, q = p — p' . The non-null diagrammatic con- 
tributions and the correspondent counterterms to 
eliminate subdivergences are presented in fig. [5J 



For the diagram of figUJi, we have 



k,i p i> - $ i> -l 



X7m 



i>-r t3 j>-k" a k 2 (l-k) 



7a" 



g ap g p ° 



(76) 



It is convenient to isolate the subdiagram, which is 
the one-loop vertex correction and contract it with 
the momentum of the external photon, q^. How- 
ever, at one- loop order, the Ward identity relating 
the vertex function and the fermionic self-energy 
is written as 

^r« = £(%)-E(V) , 

where H^{p) is the one-loop self-energy of a 
fermion with momentum p. This allows us to write 



1 



Ik i> 



:(-«0 



(-ie) Y,W(p-k)-T,W(p' -k) 



j>-f a ~k 2 ' 

and, substituting the expression for £ (1) (p), 

we 



(77) 
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have 



propagator. This amplitude is given by 



—leq 
2ie 5 b 

-2ie 5 b 



p r (2) 



-ieT 



(2) 
fib 



(2tt) 4 (p' - k) 2 k 2 
d 4 k jt-k 1 
(2tt) 4 (p - fc) 2 fc 2 



In 
In 



(p-fc) 2 
A 2 

(p' - fc) 2 
A 2 



(-ief 



-2 . 
(78) 



(-tef(-l) / TV 7cty7/3 



-ig 



0p 



(82) 

These integrals will be left unsolved, since they 

will be canceled by contributions coming from where > in the brackets, we have the expression for 



diagrams of figures [5jl and [5k, which are pre- 
sented with their correspondent counterterms. It 
is straightforward to show that 



the subdiagram. Applying the forest formula and 
using the expression for the one-loop correction to 
the photon propagator we have 



-ieT 



-2ie 5 b 
+2ie 5 b 



k [p'-k) 2 k 2 

A Qt-fl 

k (P~ k ) 2 k 2 



(p ~ kf 
A 2 

(p-fc) 2 



(2) 
fib 

A 



p 1 ' -f^-lf fc 2 



-2ie b b 



k (p'-k) 2 k 2 



In - 



A 2 

(p-fc) 2 
A 2 



--e 2 b[k a kfj - k 2 g a/3 ] 



In 



A 2 



5 

+ 3 



+2ie 5 b 1 »[lW{p)-p ti I{p)} 



and 



1 

I 2 ' 

(83) 

-2 

Contracting the amplitude with the external pho- 
(79) ton momentum q^ and using the identity $ = 
(p 1 — If) — (p 1 — ty) , we have, in terms of typical 
one and two-loop basic integrals, 



7 Mp(2) _ 
i x fie 

r-A 



-2«e 5 6 
+2ie 5 & 
= 2ie 5 6 



(/ 


-V) 


(p'- 


k) 2 k 2 




-V) 


(p- 


k) 2 k 2 




-V) 


(p- 


k) 2 k 2 



fib 



ip' - kf 
A 2 

(p' - fc) 2 
A 2 

(p' - kf 
A 2 



3 



w-m 

k k±(pi-k) 2 
k k\p-kf 



In - 



A 2 



5 
" 3 
5 
3 



- 2 
(80) 



+ -ze 4 6 7 ^ 2^/( 2 ) (p) - 2/(2) (p) + £ p j(p) 



^eV^V)-^')}- 

So, in the partial summation, 

-^(r( 2 ) + r2 + r(2)) = 

^ {4 2) (p) -pp/(p) - (^ 2 V) -p^(p'))} > 

(81) 



|ze 4 6 7 " 2K/( 2 )(p') - 27f(p') + V/(p') 



(84) 

The next amplitude is represented by the figHt 
and, following similar steps, the final result for this 
amplitude can be written as 



2ie 



-teq^ 2 



occurs the cancelation of the unsolved integrals 
corresponding to these diagrams. 

The next two-loop vertex correction is repre- 
sented by the diagram of figHb, where the sub- 
diagram is a correction for the one-loop photon 



-8ie 5 
+8ie 5 



flC 

A 



(jt-^)l-(p- fc) 



k,i k 2 l 2 (k — l) 2 (p — k) 2 (p' - I) 2 
A (jt - #)l ■ (jpf - fc) 

k,l k 2 l 2 (k~l) 2 (p'-k) 2 (p-l) 2 - 



(85) 
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To obtain the result above, it was necessary to 
make shifts and to use the fact that the integrals 
are odd in the exchange of both sign of p and p'. 
We do not need to calculate these integrals, as, as 
we will see, they will be canceled by some of the 
next contributions. 

Finally, the two last two-loop vertex corrections 
contributions are given by the figures [SJ and [5g. 
Calculating the amplitude for the first diagram, 
we have 



+8ie- 



(t-W- (p' - k) 



1 /i/CT 



-&ie 



M k 2 l 2 (k-l) 2 (p' -k) 2 {p-l) 2 

A y-wjp-k) 

u k 2 l 2 {k-l) 2 {p-k) 2 {p-l) 2 

'T 

A (p' - k) 



-Aie 



k>l k 2 l 2 (k-l) 2 (p' -k) 2 (p-l) 2 
%r{l^{p)-p,I{p)}, (86) 



where for the sake of subtracting the subdiver- 
gences we have discarded the divergent part of the 
subdiagram. For the last graph, we obtain 



(ji - fli ■ (p - k) 



k , k 2 i 2 (k-i) 2 (p-k) 2 ( P '-i) 2 



-8ie- 



(J 



(p' - k) 



u, k 2 l 2 (k — l) 2 {p' — k) 2 (p' - I) 



HgCT 



= &ie° 

' l kil k 2 l 2 (k-l) 2 (p-k) 2 (p<-l) 2 

+^e 5 br{l^(p')-P , ,np')}, (87) 

Summing the results for the graphs of figures [5fc, 
[5f and [Sg, we get 



-teq 



-Aie 



V{^ 2) (p)-Pp-T(p) 

-(i¥ ] (j?')-p'»m)} ( 88 ) 

It is now an easy task to check that 

gAT (2) =S (2)^)_ E (2)^') (gg) 

The calculation we have performed above with 
the contraction of the photon momentum with the 



photon-fcrmion vertex follows the traditional dia- 
grammatic proof of gauge invariance and, thus it 
is not surprising that it would work . Nevertheless, 
it is important to pay attention to this feature of 
Constrained Implicit Regularization: the method 
was developed to contain all the features which 
are necessary in order to preserve gauge invariance 
and dimension sensitive symmetries. The possi- 
bility of making shifts is implemented by taking 
care of all the surface terms. The maintenance of 
the specific dimension of the theory and the fact 
that the integrand is not modified guaranties the 
preservation of the vector algebra. So, the implicit 
regularization can even be a simple cutoff, since its 
bad features are under control. 

In the next section, we perform the calculation 
of the two loop QED /3-function. It is a nice test for 
the definition of counterterms in terms of our basic 
divergences. As we have seen, the CIR procedure 
has preserved the two- loop QED Ward Identities, 
which means the correct relations between QED 
renormalization constants. 



V. SPINORIAL QED TWO-LOOP 
/3-FUNCTION 

As a final checking, we are going to present the 
two-loop spinorial QED /3-function. The relation- 
ship between the bare and renormalized charge is 
given by 



Zi 



(90) 



As a consequence of the Ward Identities, Z\ = 
Z2, and it is obtained through the electromagnetic 
field counterterm, Z3. Until two loop order, this 
counterterm is written as 



-36lg(A 2 ) + H Wio9 ( A2) 



(91) 



We write eq. (|90[) in terms of the fine structure 
constant: 

a B = Z^a . (92) 

We are going to use the /3-function definition as 

a _ d(lna) _ 2A 2 da 

1 d(lnA) ~ a OX 2 ' [ ' 
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The expression for Z3 can be written in terms of 
the coupling constant as 



x (^,(^) - 36/^(A") + ^W^) j + • • • 

= 1 - c x a + c 2 a 2 + ... . (94) 

The bare coupling is independent of the scale en- 
ergy. Therefore, after performing the derivation of 
olb with respect to A 2 , the scale energy, we obtain 



da 



z; 



5A 2 3 9A 2 
The first term is given by 







(95) 



dZ^_ 
d\ 2 



da dci da 



After some algebra, we obtain 



(96) 



-2A 2 



[1 - 2aci +3a 2 4 H ] 

■ 8C2 



9\2 # c l 

2A a— 2 



d\ 2 



(97) 



We can write 
= 2A 2 a|^ - 2AV 



dc 2 
9A 2 



2c w 



The derivatives of c\ and C2 are given by 



da _ d 
5A 2 = 9A 2 



167Ti/, ofl (A 2 ) 



1 



3ttA 2 



(98) 
(99) 



and 
<9c 2 



If? ^ 

dA 2 ~ 3 n 1 12^ 2 A 2 



which leads us to 

dc 2 _ 9ci_ 
5A 2 X 9A 2 



WA 2 ) 



2tt 4 A 2 (16) 2 

(100) 



-16tt 2 



2tt 4 A 2 (16) 2 



9A 2 } 



4ttA 2 



(101) 



Finally, we can write the result for the [3 function 
obtained up to two loop order: 



2a 1 /a\ 2 
~3n + 2 



(102) 



which is in agreement with the known result, since 
the (3 function is independent of the renormaliza- 
tion scheme up to two loop order. 



VI. CONCLUSION 

In this paper we extended the procedure of Con- 
strained Implicit Regularization to be applied in 
abclian gauge theories beyond one-loop order. For 
the sake of simplicity, we have considered the mass- 
less case. The massive case does not present any 
new features, apart from the fact that the calcu- 
lation is a little more involved. The whole process 
is based on the elimination of symmetry violating 
terms, which are momentum space surface terms. 
This is equivalent to allow shifts in the momen- 
tum of integration. This is one of the principal 
requirements in the diagrammatic proof of gauge 
invariance and is the reason why Dimensional Reg- 
ularization preserves this symmetry. 

In the case of anomalies, the elimination of 
the surface terms by means of symmetry restor- 
ing counterterms does not work. The reason is 
that the counterterm necessary to restore, for in- 
stance, the vector Ward identity would cause the 
violation of the axial Ward identity and vice- versa. 
Actually, this occurs because, diagrammatically, 
the anomaly manifests itself as a violation of mo- 
mentum routing invariance. For methods like Di- 
mensional Regularization, which eliminate surface 
terms automatically, the ambiguity is manifested 
by the form the 75 matrix is inserted into the trace. 

For the non-abelian case, CIR was successfully 
applied in the renormalization on QCD to one-loop 
order 12]. We argue here that it is also adequate 
for higher order calculations. Our argument is 
based on the fact that the 't Hooft identity [27| re- 
quires a regularization method which allows shifts 
in the integration momenta. This is well discussed 
in the references [28[ and [29j. Implicit Regular- 
ization contains all the ingredients to preserve the 
generalized Ward-Slavnov- Taylor identities and we 
expect it would succeed in higher order calcula- 
tions also in the non-abelian case as well. That is 
because after group theoretical factors have been 
evaluated, we are left with the same basic integrals, 
just as in the abelian case. 

Actually, the main utility of Implicit Regular- 
ization is its application in dimension specific theo- 
ries such as supersymmetric theories. As discussed 
in the introduction, from a more phenomcnolog- 
ical perspective it is expected that the Minimal 
Supersymmetric Standard Model (MSSM) can be 
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probed at LHC with an accuracy at the percent 
level in measurements of electroweak precision ob- 
servables Q. The theoretical evaluation of such 
observables must be performed at least to two 
loop order both in the SM and MSSM so that 
the evidence of new physics becomes unraveled. 
Clearly, an invariant regularization and renorma- 
lization procedure is important, both to eliminate 
inconsistencies and to reduce the number of ad- 
justable parameters in the calculation. 



VII. APPENDIX 

In this appendix we display the results of the 
basic integrals up to two loop order and illustrate 
the procedure with the evaluation of one of them, 
including the finite part, which is similar for the 
others integrals. We begin by the explicit calcula- 
tion of the two-loop integral, 

V' 1 ' = f A ^— t1>. |-^|. (103) 



J , ( Jf_ 

k k*{k- P y n \ a 2 

in which we will insert a fictitious mass m that will 
be set to zero in the end of the calculation. The 
use of the identity JT]) gives us 



j(2) 



A 



1 



hi 



k (fc 2 — TO 2 ) 2 

2 - 2p ■ k 



m 



(fc 2 — m 2 ) 2 [(p — fc) 2 — m 2 



In 



A 2 



4> 2 w (2) - 



(104) 

It is clear from the equation above the reason for 
the introduction of the fictitious mass m. Although 
the integral by itself is infrared finite, when the sep- 
aration by means of the relation |T]) is performed, 
we are left with two infrared divergent parts. The 
scale relation, 



:- 



will make the connection between these parts in 
order to make the limit m 2 — > possible. We show 
a simple trick to calculate the finite part. First, we 
make use of the standard limit 

a e - 1 



to write 



In a = lim ■ 



/(^lim-Vi, 

1=0 



(105) 



(106) 



with 

A, = 



(-1) 



i-/ 



p 2 — 2p ■ k 



1 (-A 2 ) k J k (fc 2 - m 2 ) 2 - u [{p - fc) 2 - m 2 ] ' 

(107) 

The A[ terms can be calculated by means of Feyn- 
man parametrization. As it is easy to note, only 
the coefficient of e in the series for e — > of A\ will 
contribute. So, we get, for small m, 

,2 s 



p) = -b 








-M- 


p Z 




A 2 



m: 
A 2 " 



m 
A 2 " 



o 

p 
A 2 



(108) 



When the divergent and finite parts are added to- 
gether and the scale relation is used, the depen- 
dence on the fictitious mass m disappears and we 
get 

^>=^) + 4X4) +!»(-£ 

(109) 

The results of the other integrals, as defined in 
the text, are the following: 



A 2 



I = Ii og (X 2 )-bln 



26; 



2b 



(110) 



(111) 



V = \vi>,Vv \liog(X 2 ) - bin ) + ^-b 



1 2 J 1 r l\2\ , b i„ I P" 



+ -p Z 9^\—^o g {^) + -^[- x2 



4 2) = y{4 2 j(A 2 ) + ^ 03 (A 2 ) 



(112) 



(113) 



^ [ -~ 2 



1 / p 2 

3 ln r^ 



/ '"'' / "' 1 -45(A 2 ) + ^ 09 (A 2 )+ 



12 

1 ^ 2 
2 ln 



(U4) 



16 



J(2) = 

+3fe] 



o2 r 



I log {\ 2 )-b\n 



Hi 

A 2 



(115) 



For the finite part, we adopt the following proce- 
dure. We write 



B^ v = Cip^pv + C 2 p 2 gpLv 



(118) 



Another class of integrals we have to deal to two 
loop order are the overlapped ones. In this case, 
there is a symmetry between the internal momen- 
tum k and I in the denominator, with three factors 
involving each one. This type of integrals com- 
monly appears in the solution of overlapped dia- 
grams. So it is interesting to maintain all the sub- 
divergences to be eliminated at the end, when the 
counterterms are added. We give some details of 
the solution of one of them: 



jOl = f 

Jk 



k^k v 



k,l k 2 l 2 {k-l) 2 (p-l) 2 (p-kf 

k^i ki/ 



Ik, 



k,l fc 4 / 2 (fc-/) 2 (p-0 2 

(p 2 - 2p ■ fykpkv 



k 4 l 2 (k~l) 2 (p-l) 2 {p-k) 2 



-A — R 



(116) 



and by contracting the integral with p^Pv and 
with g^ v , we solve a system of two equations to ob- 
tain C\ and C2. This procedure is not allowed for 
divergent integrals in the context of Implicit Reg- 
ularization, since the contraction with g^ u would 
make us loose control of the surface terms that 
must be eliminated. The final result for I®J yields 



C = 9 ~f {^(A 2 ) - <(A 2 ) + 9 -bI log (X 2 ) 
W-g)WA 2 ) + ^ln 2 (-£ 



7 



b 2 In 



p_ 

A 2 



^ 31 
+ — 



b 2 - V -I° 



(119) 



The second integral is finite. The first one is easily 
solved: 



A f A 1 f A k^ky 
» v J l l 2 (p-l) 2 )J k k\k- 

r A i_ 

~Ji W z 



g^v 
i) 2 ) \ T 



-i) 2 

I log (X 2 )-b\n(- 



+ m + b - l f 



= ^(WA 2 ) + ^)/-&^/( 2 ) 
b 2 p^p v 



+ 



4 p 2 



A 2 



(117) 



For the other divergent integral of this type we 
needed, we have: 



7-02 



l ±b 2 -b 2 *- 

12 36 



+ 



g»v \ ^Ilog(X Z ) - T ln ( ~ 



12 



(120) 
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